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Abstract 

The determination of twist-4 corrections to the structure functions of polarized 
e(n)N scattering by QCD sum rules is reviewed and critically analyzed. It is found 
that in the case of the Bjorken sum rule the twist-4 correction is small at 
Q 2 > 5 GeV 2 . However, the accuracy of the today experimental data is insufficient 
to reliably determine a s from the Bjorken sum rule. For the singlet sum rule - p + n 
- the QCD sum rule gives only the order of magnitude of twist-4 correction. At low 
and intermediate Q 2 the model is presented which realizes a smooth connection of 
the Gerasimov-Drell-Hearn sum rules at Q 2 = with the sum rules for T Ptn (Q 2 ) at 
high Q 2 . The model is in a good agreement with the experiment. 



1 



1. Introduction 

In the last few years there is a strong interest to the problem of nucleon spin structure: 
how nucleon spin is distributed among its constituents - quarks and gluons. New experi- 
mental data continuously appear and precision increases (for the recent data see [1], [2]). 
One of the most important item of the information comes from the measurements of the 
first moment of the spin-dependent nucleon structure functions g\ (x) which determine the 
parts of nucleon spin carried by u, d and s quarks and gluons. The accuracy of the data 
is now of a sort that the account of twist-4 terms is of importance when comparing the 
data with the Bjorken and Ellis- Jaffe sum rules at high Q 2 . On the other side, at low and 
intermediate Q 2 a smooth connection of the sum rules for the first moments of gi(x, Q 2 ) 
with the Gerasimov-Drell-Hearn (GDH) sum rules [3,4] is theoretically expected. This 
connection can be realized through nonperturbative Q 2 -dependence only. In my talk I 
discuss such nonperturbative <5 2 -dependence of the sum rules (see also [5]). 

Below I will consider only the first moment of the structure function g±(x, Q 2 ) 

i 

r p ,„(Q 2 ) = y dx gi; P , n (x,Q 2 ) (1) 



The presentation of the material is divided into two parts. The first part deals with 
the case of high Q 2 . I discuss the determination of twist-4 contributions to r Pi „ by QCD 
sum rules and with the account of twist-4 corrections and of the uncertainties in their 
values compare the theory with the experiment. In the second part the case of low and 
intermediate Q 2 < lGeV 2 is considered in the framework of the model, which realizes the 
smooth connection of GDH sum rule at Q 2 = with the asymptotic form of Y p ^ n {Q 2 ) at 
high Q 2 . 

2. High Q 2 

At high Q 2 with the account of twist-4 contributions Y p ^ n {Q 2 ) have the form 

r P AQ 2 ) = r a P %(Q 2 ) + K w n \Q 2 ) (2) 

+ i[l - l -a - 0.55a 2 - 4.45a 3 ]s| - ^-a s (Q 2 )Ag(Q 2 ) (3) 

In eq.(3) a = a s (Q 2 )/-K, g A is the /3-decay axial coupling constant, g A = 1.260 ± 0.002 [6] 

g A = Au- Ad a 8 = Au + Ad - 2As E = Au + Ad + As. (5) 

Au, Ad, As, Ag are parts of the nucleon spin projections carried by u, d, s quarks and 
gluons: 
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dx (6) 

where q+(x), q~(x) are quark distributions with spin projection parallel (antiparallel) to 
nucleon spin and a similar definition takes place for Ag. The coefficients of perturbative 
series were calculated in [7-10], the numerical values in (3) correspond to the number of 
flavours Nf = 3, the coefficient c was estimated in [11], c ~ 130. In the MS renormal- 
ization scheme chosen in [7-10] a 8 and E are Q 2 -independent. In the assumption of the 
exact SU(3) flavour symmetry of the octet axial current matrix elements over baryon 
octet states a 8 = 3F - D = 0.59 ± 0.02 [12]. 

Strictly speaking, in (3) the separation of terms proportional to S and Ag is arbitrary, 
since the operator product expansion (OPE) has only one singlet in flavour twist- 2 op- 
erator for the first moment of the polarized structure function - the operator of singlet 

axial current j^(x) = J2 Qi{ x )lfil5q, q — u,d,s. The separation of terms proportional 

q 

to £ and Ag is outside the framework of OPE and depends on the infrared cut-off. The 
expression used in (3) is based on the physical assumption that the virtualities p 2 of gluons 
in the nucleon are much larger than light quark mass squares, \p 2 \ ^> m 2 [13] and that 
the infrared cut-off is chosen in a way providing the standard form of axial anomaly [14]. 

Twist-4 corrections to r Pi „ were calculated by Balitsky, Braun and Koleshichenko 
(BBK) [15] using the QCD sum rule method. 

BBK calculations were critically analyzed in [5], where it was shown that there are 
few possible uncertainties in these calculations: 1) the main contribution to QCD sum 
rules comes from the last accounted term in OPE - the operator of dimension 8; 2) there 
is a large background term and a much stronger influence of the continuum threshold 
comparing with usual QCD sum rules; 3) in the singlet case, when determining the induced 
by external field vacuum condensates, the corresponding sum rule was saturated by rj- 
meson, what is wrong. The next order term - the contribution of the dimension 10 
operator in the BBK sum rules was estimated by Oganesian [16]. The account of the 
dimension-10 contribution to the BBK sum rules and estimation of other uncertainties 
results in (see [5]): 

b p _ n = -0.006 ± 0.012 GeV 2 (7) 

b p+n = -0.035(±100%) GeV 2 (8) 

As is seen from (7), in the nonsinglet case the twist-4 correction is small (< 2% at 
Q 2 > 5GeV 2 ) even with the account of the error. In the singlet case the situation is much 
worse: the estimate (8) may be considered only as correct by the order of magnitude. 

I turn now to comparison of the theory with the recent experimental data. In Table 1 
the recent data obtained by SMC [1] and E 154(SLAC) [2] groups are presented. 



A, = /, 



q+(x) - q-(x) 
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Table 1 











a s (5GeV z ) 


SMC 


0.132 ±0.017 


-0.048 ± 0.022 


0.181 ±0.035 


o.270^:iS 


combined 


0.142 ±0.011 


-0.061 ±0.016 


0.202 ±0.022 


0.1161^ 


E 154(SLAC) 


0.112 ±0.014 


-0.056 ±0.008 


0.168 ±0.012 


0.339 ™ 


EJ/Bj sum rules 


0.168 ±0.005 


-0.013 ±0.005 


0.181 ±0.002 


0.276 



In the second line of Table 1 the results of the performed by SMC [1] combined 
analysis of SMC [1], SLAC-E80/130 [17], EMC [18] and SLAC-E143 [19] data are given. 
The data presented in Table 1 refer to Q 2 = 5GeV 2 . In all measurements each range 
of x corresponds to each own mean Q 2 . Therefore, in order to obtain gi(x,Q 2 ) at fixed 
Q 2 refs. [1,2] use the following procedure. At some reference scale Qq (Qq = lGeV 2 
in [1] and Ql = 0.34:GeV 2 in [2]) quark and gluon distribution were parametrized as 
functions of x. (The number of the parameters was 12 in [1] and 8 in [2]). Then NLO 
evolution equations were solved and the values of the parameters were determined from 
the best fit at all data points. The numerical values presented in Table 1 correspond to 
MS regularization scheme, statistical, systematical, as well as theoretical errors arising 
from uncertainty of a s in the evolution equations, are added in quadratures. In the last 
line of Table 1 the Ellis- Jaffe (EJ) and Bjorken (Bj) sum rules prediction for T p , T n and 
T p — T n , correspondingly are given. The EJ sum rule prediction was calculated according 
to (3), where As = , i.e., £ = as = 0.59 was put and the last-gluonic term in (3) was 
omitted. The twist-4 contribution was accounted in the Bj sum rule and included into 
the error in the EJ sum rule. The a s value in the EJ and Bj sum rules calculation was 
chosen as a s (5GeV 2 ) = 0.276, corresponding to a s (M z ) = 0.117 and A^L = 3Q0MeV (in 
two loops). As is clear from Table 1, the data, especially for T n , contradict the EJ sum 
rule. In the last column, the values of a s determined from the Bj sum rule are given with 
the account of twist-4 corrections. 

The experimental data on T p presented in Table 1 are not in a good agreement. 
Particularly, the value of T p given by E154 Collaboration seems to be low: it does not 
agree with the old data presented by SMC [20] (T p = 0.136 ± 0.015) and E143 [19] 
(r p = 0.127±0.011). Even more strong discrepancy is seen in the values of a s , determined 
from the Bj sum rules. The value which follows from the combined analysis is unacceptably 
low: the central point corresponds to A^± = 15MeV\ On the other side, the value, 
determined from the E154 data seems to be high, the corresponding a s (M z ) = 0.126 ± 
0.009. Therefore, I come to a conclusion that at the present level of experimental accuracy 
a s cannot be reliably determined from the Bj sum rule in polarized scattering. 

Table 2 shows the values of E - the total nucleon spin projection carried by u, d and 
s-quarks found from Y p and T n presented in Table 1 using eq.(3). (It was put Qa = 
1.260, as = 0.59, the term, proportional to Ag is included into E.). 
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Table 2: The values of £ 





From T p 


From T n 




At a s (5GeV 2 ) = 
= 0.276 


At a s (GeV 2 ) 
given in Table 1 


At a s (5GeV 2 ) = 
= 0.276 


At a s (5GeV 2 ) 
given in Table 1 


SMC 


0.256 


0.254 


0.264 


0.266 


Comb. 


0.350 


0.250 


0.145 


0.225 


E154 


0.070(0.14; 0.25) 


0.133(0.20; 0.30) 


0.19(0.25; 0.14) 


0.144 (0.205; 0.10) 



In their fitting procedure [2] E154 Collaboration used the values a$ = 0.30 and qa = 
1.09. The values of £ obtained from T p and T n given by E154 at a 8 = 0.30, = 1.26 
and as = 0.30, 5u = 1.09 are presented in parenthesis. The value as = 0.30 corresponds 
to a strong violation of SU(3) flavour symmetry and is unplausible; qa = 1.09 means a 
bad violation of isospin and is unacceptable. As seen from Table 1, E is seriously affected 
by these assumptions. The values of £ found from T p and T n using SMC and combined 
analysis data agree with each other onlyif one takes for a s (5GeV 2 ) the values given in 
Table 1 (ct s = 0.116 for combined data), what is unacceptable. The twist-4 corrections 
were not accounted in E in Table 2: their account, using eqs.(7), (8), results in increasing 
of £ by 0.04 if determined from T p and by 0.03 if determined from T n . 

To conclude, one may say, that the most probable value of E is E m 0.3 with an 
uncertain error. The contribution of gluons may be estimated as Ag(lGeV 2 ) ~ 0.3 (see 
[5], [21], [22]). Then Ag{hGeV 2 ) ~ 0.6 and the account of gluonic term in eq.(3) results 
in increasing of E by 0.06. At E = 0.3 we have Au = 0.83, Ad = -0.43, As = -0.1. 

3. Low and Intermediate Q 2 

The problem of a smooth connection of the Gerasimov-Drell-Hearn (GDH) sum rules 
[3,4] which holds at Q 2 = 0, and the sum rules at high Q 2 attracts attention in the last 
years [5, 23-25]. 

In order to connect the GDH sum rule with T p ^ n (Q 2 ) consider the integrals [26] 

WQ 2 )= / ^G up , n (v,Q 2 ) (9) 

Q 2 /2 

Changing the integration variable v to x, (9) can be also identically written as 
At Q 2 = the GDH sum rule takes place 



W) = —^l = -0.8035; J n (0) = --/£== 0.9149; I p (0) - J n (0) = 0.1114 (11) 

where k p and proton and neutron anomalous magnetic moments. 

The schematic Q 2 dependence of I P (Q 2 ), I n {Q 2 ) and I P (Q 2 ) — I n {Q 2 ) is plotted in 
Fig.l. The case of I p {Q 2 ) is especially interesting: I p (Q 2 ) is positive, small and decreasing 
at Q 2 > 3GeV 2 and negative and relatively large in absolute value at Q 2 = 0. With 
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I n (Q 2 ) the situation is similar. All this indicates large nonperturbative effects in I(Q 2 ) 
at Q 2 < lGeV 2 . 

In [23] the model was suggested, which describes I{Q 2 ) (and T(Q 2 )) at low and inter- 
mediate Q 2 , where GDH sum rules and the behaviour of I(Q 2 ) at large Q 2 where fullfilled. 
The model had been improved in [24]. (Another model with the same goal was suggested 
by Soffer and Teryaev [25]). 

Since it is known, that at small Q 2 the contribution of resonances to I{Q 2 ) is of 
importance, it is convenient to represent I{Q 2 ) as a sum of two terms 

I(Q 2 ) = I res (Q 2 ) + I\Q 2 ), (12) 

where I res (Q 2 ) is the contribution of baryonic resonances. I res (Q 2 ) can be calculated from 
the data on electroproduction of resonances. Such calculation was done with the account 
of resonances up to the mass W = 1.8GeV [27]. 

In order to construct the model for nonresonant part I'{Q 2 ) consider the analytical 
properties of I(Q 2 ) in q 2 . As is clear from (9), (10), I(Q 2 ) is the moment of the structure 
function, i.e. it is a vertex function with two legs, corresponding to ingoing and outgoing 
photons and one leg with zero momentum. The most convenient way to study of analyt- 
ical properties of I{q 2 ) is to consider a more general vertex function I{q 2 -,q 2 \p 2 ), where 
the momenta of the photons are different, and go to the limit p — > 0, q\ — > q\ = q 2 . 
I{q 2 ,q2',P 2 ) can be represented by the double dispersion relation: 

I(q 2 ) = n lim l{qlq 2 2] p 2 ) = [( dsj ds, - P^lfl ■ 



«^=^o > {J ^ (ai-gf)( S2 -d) 

f^4ds + P(qi) f^4ds) (13) 

The last two terms in (13) are the substruction terms in the double dispersion relation, 
P(q 2 ) is the polynomial. According to (10), I{q 2 ) decreases at | q 2 |— > oo, P(q 2 ) = 
Const and the constant subtraction term in (13) is absent. We are interesting in I{Q 2 ) 
dependence in the domain Q 2 < lGeV 2 . Since after performed subtraction, the integrals 
in (13) are well converging, one may assume, that at Q 2 < 2— 3GeV 2 the main contribution 
comes from vector meson intermediate staties, so the general form of I'(Q 2 ) is 

where A and B are constants, p, is p (or uj) mass. The constants A and B are determined 
from GDH sum rules at Q 2 = and from the requirement that at high Q 2 ^> p 2 takes 
place the relation 

i{Q 2 ) » i\q 2 ) » ^r-(g 2 ), (is) 

where T as (Q 2 ) is given by (3). (I res (Q 2 ) fastly decreases with Q 2 and is very small above 
Q 2 = 3GeV 2 ). These conditions are sufficient to determine in unique way the constant A 
and B in (14). For I'{Q 2 ) it follows: 
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I'(Q 2 ) = 2m 2 T as (Q 2 J 



1 



(16) 



Q2 +/i 2 (g2 +/i 2)2 



1 1 



K 2 + r es (o) , 



(17) 



c= 1 + 



2m 2 T as (Q 2 J [4 



where /; es (0) = -1.03, /; es (0) = -0.83 [24]. 

The model and eq.16 cannot be used at high Q 2 > 5GeV 2 : one cannot believe, that 
at such Q 2 the saturation of the dispersion relation (13) by the lowest vector meson is 
a good approximation. For this reason there is no matching of (16) with QCD sum rule 
calculations of twist-4 terms. (Formally, from (16) it would follow 6 p _ n rs —0.15, b p+n ps 
—0.07). It is not certain, what value of the matching point Q 2 m should be chosen in 
(16). This results in 10% uncertainty in the theoretical predictions. Fig.'s 2,3 shows the 
predictions of the model in comparison with recent SLAC data [28] , obtained at low Q 2 = 
0.5 and 1.2GeV 2 as well as SMC and SLAC data at higher Q 2 . The chosen parameters are 
T p {Q 2 m ) = 0.142, T n {Q 2 m ) = -0.061, corresponding to c p = 0.458, c n = 0.527 in (16), (17). 
The agreement with the data, particularly at low Q 2 , is very good. The change of the 
parameters only weakly influences Y p ^ n (Q 2 ) at low Q 2 . 

This work was supported in part by the Russian Foundation of Fundamental Re- 
search Grant 97-02-16131, CRDF Grant RP2-132 and Schweizerische Nationalfond Grant 
7SUPJ048716. 
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Figure Captions 

Fig. 1 The Q 2 -dependence of integrals I P (Q 2 ) , I n (Q 2 ) , I p (Q 2 ) - I n (Q 2 ). 
The vertical axis is broken at negative values. 

Fig. 2 The Q 2 -dependence of T p = T' p + IT es (solid line), described by 
eqs.(12,16,17). T p es ' (dotted) and T' p (dashed) are the resonance 
and nonresonance parts. The experimental points are: the dots 
from E143 [28], the square - from E143 (SLAC) [19], the cross - 
SMC-SLAC combined data [1], the triangle from SMC [1]. 

Fig. 3 The same as in Fig.2,but for neutron. The experimental points are: 
the dots from E143 (SLAC) measurements on deuteron [28], the 
square at Q 2 = 2GeV 2 is the E142(SLAC) [29] data from measure- 
ments on polarized 3 He, the square at Q 2 = 3GeV 2 is E143(SLAC) 
[30] deuteron data, the cross is SMC-SLAC combined data [1], the 
triangle is SMC deuteron data [1]. 
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Figure 1: 
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Figure 2: 
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Figure 3: 
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